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ABSTRACT 

We present results of variability search in the field of the young open cluster NGC 1502. Eight 
variable stars were discovered. Of six other stars in the observed field that were suspected for vari- 
ability, we confirm variability of two, including one (3 Cep star, NGC 1502-26. The remaining four 
suspects were found to be constant in our photometry. In addition, U BVIq photometry of the well- 
known massive eclipsing binary SZ Cam was obtained. 

The new variable stars include: two eclipsing binaries of which one is a relatively bright detached 
system with an EA-type light curve, an a 2 CVn-type variable, an SPB candidate, a field RR Lyr star 
and three other variables showing variability of unknown origin. The variability of two of them is 
probably related to their emission in H a , which has been measured by means of the a index obtained 
for 57 stars brighter than V ~ 16 mag in the central part of the observed field. Four other non- variable 
stars with emission in Ha were also found. 

Additionally, we provide VIq photometry for stars down to V = 17 mag and UB photometry 
for about 50 brightest stars in the observed field. We also show that the 10-Myr isochrone fits very 
well the observed color-magnitude diagram if a distance of 1 kpc and mean reddening, E(V — Iq) = 
0.9 mag, are adopted. 

Key words: stars: early-type - stars: emission-line, Be — open clusters and associations: individ- 
ual: NGC 1502 



1. Introduction 

This is the seventh paper in the series in which we present results of the vari- 
ability search for B-type variables (mainly [3 Cep, SPB and massive binaries) in 
the young open clusters of the Northern Hemisphere. The results for previously 
studied clusters were briefly summarized by Pigulski et al. (2002). In addition to 
the discovery or verification of variability of about 70 stars, mostly members of the 
investigated clusters, it appeared that the incidence of [3 Cep stars in the northern 
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clusters is lower than in the southern ones. This was interpreted in terms of the 
metallicity gradient in the Galaxy. One of the main motivations for the search is 
also a selection of clusters containing large numbers of pulsators suitable for aster- 
oseismology. An example is the previous paper of the series (Kolaczkowski et al. 
2004), in which we investigated the open cluster NGC 6910 discovering four [3 Cep 
stars. This prompted us to choose NGC 6910 as a target of a two-year photomet- 
ric and spectroscopic campaign. The data from the campaign are in the course of 
reduction, but preliminary analysis already increased the number of [3 Cep stars in 
this cluster to seven (Pigulski et al. 2007) making it an excellent target for seismic 
modelling and one of very few which are rich in [3 Cep stars. 

In the present paper, we report the results of the variability search in the field 
of another young open cluster, NGC 1502. For the first time, the data from the new 
CCD camera, with a much larger field of view, are included. 

2. The Cluster 

NGC 1502 (C 0403 +622, OC1383) is a loose young open cluster of Trumpler 
type II3p in Camelopardalis. Almost all studies locate the cluster near the outer 
border of the Orion Spur at a distance of 0.7-1.0 kpc (Sanford 1949, Johnson et 
al. 1961, Purgathofer 1961, Nicolet 1981, Reimann and Pfau 1987, Kharchenko et 
al. 2005, Paunzen et al. 2005), and only Tapia et al. (1991) derive a significantly 
larger distance of 1.5 kpc. Because the main sequence of the cluster does not have 
a well-defined turnoff, the estimated age of NGC 1502 has a large spread, ranging 
from 5 to 15 Myr (Lindoff 1968, Harris 1976, Mermilliod 1981, Reimann and Pfau 
1987, Tapia et al. 1991, Kharchenko et al. 2005). Due to the similarity of ages 
and distances, a nearby 6° apart) Cam OBI association is frequently related to 
NGC 1502 (de Zeeuw et al. 1999, Lyder 2001, Straizys and Laugalys 2007) and it 
is possible that both objects have a common origin. It was also suggested (Blaauw 
1954, 1961, 1964; Stone 1991) that the 09.5 Iae supergiant a Cam is a runaway 
star from NGC 1502. The proper motions of stars in the field of the cluster were 
considered in many studies devoted to the kinematics of the Galactic system of open 
clusters (see e.g., Rastorguev et al. 1999), but some, like those of de Vegt (1966) or 
Dias et al. (2001), may also be used to investigate the cluster membership. 

The mean E(B — V) color-excess of members of NGC 1502 amounts to about 
0.75 mag (e.g., Reimann and Pfau 1987) but seems to vary in a range of the order 
of 0.2 mag across the cluster (Janes and Adler 1982, Yadav and Sagar 2001). How- 
ever, as found by Tapia et al. (1991) and confirmed by Pandey et al. (2003) and 
Weitenbeck et al. (2008), the total-to-selective absorpion ratio Ry = Ay /E(B — V) 
is for this region of the sky significantly smaller than the mean Galactic value 
of 3.1 ± 0.1. The three studies provide Ry between 2.4 and 2.7. The study of 
Weitenbeck et al. (2008) was based on polarimetric observations. They found that 
polarization for the cluster stars slightly exceeds 6 %, in accordance with earlier 
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measurements (Hall and Mikesell 1950, Dombrovskii and Hagen-Thorn 1964). 

The field of NGC 1502 has been a subject of many photometric studies. The 
first two-color photographic photometry was published by Zug (1933). The UBV 
photometry, both photographic and photoelectric, was carried out by Hoag et al. 
(1961), Purgathofer (1964), Dombrovskii and Hagen-Thorn (1964) and Reimann 
and Pfau (1987). Stromgren uvby photometry was obtained by Reimann and Pfau 
(1987), Tapia et al. (1991), Delgado et al. (1992) and Crawford (1994). The three 
latter papers contain also fj photometry for some cluster members. The near- 
infrared JHK photometry was provided by Tapia et al. (1991) and observations in 
the narrow-band Aa photometric system were published by Paunzen et al. (2005). 
The most thorough spectroscopic survey of the cluster was done by Huang and Gies 
(2006ab). They published MK spectral types, projected rotational velocities, ra- 
dial velocities for two epochs, effective temperatures, surface gravities and helium 
abundances for 18 stars in the field of NGC 1502. Some other authors published 
MK spectral types and/or radial velocities for several bright cluster members (e.g., 
Hoag and Applequist 1965). The mean radial velocity of the cluster amounts to 
-21 ± 6kms"' (Kharchenko et al. 2005). 

There are several numbering systems used for stars in the field of NGC 1502 
(Zug 1933, Hoag et al. 1961, Purgathofer 1964, Dombrovskii and Hagen-Thorn 
1964, Tapia et al. 1991). Unless noted, we will follow the numbering provided by 



The most extensively studied star in NGC 1502 is the eclipsing binary SZ Cam, 
the north-western component of the visual binary ADS 2984, the bright pair in the 
center of the cluster (Fig.Q}. SZ Cam = HD 25638 = ADS 2984B is itself a quadru- 
ple hierarchical system consisting of two spectroscopic binaries of which one is 
eclipsing. The system will be presented in more detail in Section 4.3. The other 
component of the pair, ADS 2948A = HD 25639, is also a spectroscopic binary. 
Several other stars in the field of NGC 1502 were claimed to be variable or were 
suspected for variability. Guthnick and Prager (1930) noted photometric variabil- 
ity of HD 25639 = NGC 1502-1 with an amplitude of 0.05-0.06 mag. Next, two 
stars measured by Purgathofer (1964), NGC 1502-44 (designated F in his paper) 
and 5, were labeled by him as 'var?'. Hill (1967) observed photoelectrically six 
stars in NGC 1502 and found variability with the B amplitude of 18 mmag and 
a period of 0.190 d or 0.209 d in the combined light of one of them, namely a 
close visual pair, NGC 1502-37 and 38. He concluded that one of the stars in the 
pair is a p Cep-type candidate. Subsequently, from Stromgren photometry car- 
ried out on several nights, Delgado et al. (1992) found three variable candidates in 
NGC 1502 among 22 stars checked for variability. These were NGC 1502-1, 26 and 
61. NGC 1502-26, showing frequencies of 10.5 d -1 and 15.7 d _1 , appeared to be 
a very good [3 Cep candidate. Another [3 Cep candidate was NGC 1502-1. For this 
star, Delgado et al. (1992) found periodic variations with a frequency of 5.2 d _1 

' The WEBDA database is available at |http : //www.univie ■ ac ■ at/webda/webda .html| 
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superimposed on a variability on a longer time scale, attributed to its binarity. The 
third star, NGC 1502-61, they suspected to be variable on a long time scale. 
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Fig. 1. A 20' x 20' fragment of the POSS-II infrared plate centered on NGC 1502. Solid lines delimit 
the fields observed with the two CCD cameras we used, see text for explanation. North is up, East, to 
the left. 

The CCD photometry of NGC 1502 we present here allows us to verify variabil- 
ity of all suspects and candidates mentioned above and carry out a comprehensive 
variability study of stars in the cluster field. The results are discussed in Section 4; 
some preliminary findings were already published by Steslicki (2006). 

3. Observations and Reductions 

All observations of NGC 1502 were carried out at Bialkow, a station of Uni- 
versity of Wroclaw, with a 60-cm Cassegrain telescope. During the first observing 
run, we acquired observations on 22 nights between August 22, 2004 and April 
20, 2005. We used a 576 x 384 pixels Star I CCD camera, the same which was 
used to obtain the data for other clusters that were the subjects of previous papers 
of the series. The camera covered a rectangle of only 6'x 4'. For this reason, we 
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observed two overlapping fields to obtain photometry of a larger number of cluster 
members. These fields are the two overlapping rectangles shown in Fig. [T] The 
background image comes from the POSS-ltj infrared plate centered on NGC 1502. 

Starting from the late summer of 2005, the observations of NGC 1502 were re- 
sumed with the new, Andor DW432 B V CCD camera having a 1250 x 1 152 pixels 
back-illuminated chip. The use of the new camera resulted in a larger field of view, 
comprising about 13' x 12' (large rectangle in Fig. Q]). These observations were 
made on three nights, between September 8 and October 7, 2005. In addition, the 
cluster was observed on 12 nights between January 3 and March 30, 2008. In total, 
we gathered almost 14 000 frames, which resulted in 162 hours of time-series data 
(90 hours with old and 72 hours with the new camera). Observations were made 
primarily in Johnson V and Cousins Iq -bands. In 2008, some frames through John- 
son U and B filters were also taken. Finally, like for most clusters we observed so 
far, Ha images were secured for the purpose of finding stars with emission in this 
spectral line (see, e.g., Kolaczkowski et al. 2004). 

Depending on the camera and the sky conditions, the exposure times for time- 
series data were equal to 20-50 s for V and 20-30 s for Iq ■ In order to increase 
dynamical range of the photometry, shorter exposures were also taken from time to 
time. This was done mainly in order to get better photometry for the two brightest 
stars in the field, the visual double ADS 2984. All frames were calibrated in a 
standard way and reduced by means of the DAOPHOT II package (Stetson 1987). 
Both profile and aperture photometries were obtained, the latter by means of the 
DAOGROW program (Stetson 1990). We used both types of photometry: usually 
for bright unsaturated stars the aperture photometry appeared to have smaller scatter 
than the profile photometry, while for the faint stars the opposite was true. 

4. Variable Stars 

As noted above, Steslicki (2006) already presented the results of the analy- 
sis of almost all frames of NGC 1502 acquired with the old camera. These re- 
sults included the new data for SZ Cam = NGC 1502-2, confirming variability of 
NGC 1502-26, the [3 Cep suspect, and a discovery of two other variable stars, 
NGC 1502-7 and 43. The former star was found to be eclipsing with an orbital 
period of 1.9839 d, while the variability of the latter was best represented by a 
sinusoid with a period of 0.41885 d. 

The observations obtained with the new CCD camera not only cover much 
larger area of the cluster and its surroundings but also have a better quality. For 
this reason, in the search for variability we used primarily the data from the new 
camera. However, a search for variable stars was also carried out using the data 

^The Second Palomar Observatory Sky Survey (POSS-II) was made by the California Institute of 
Technology with funds from the National Science Foundation, the National Geographic Society, the 
Sloan Foundation, the Samuel Oschin Foundation, and the Eastman Kodak Corporation. 
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obtained with the old camera in order to be sure that no transient changes or eclipses 
were missed. In the final analysis, data from both cameras were combined for stars 
found to be variable. 
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Fig. 2. Schematic map of the observed field with positions of all stars detected in the Iq band. 
Variable stars are marked and labeled (see Table 1 for designations). North is approximately up, East, 
to the left. 



Comparison stars were chosen in a trial-and-error procedure out of the bright 
unsaturated stars. The differential magnitudes were analyzed by means of a Fourier 
periodogram routine in the range between and 40 d ~ 1 . Then, the periodograms, 
the light curves and phase diagrams were inspected by eye and the final decision as 
to the variability was taken. Among 1034 stars detected in the field of view of the 
new camera eleven showed variability. They are indicated in the schematic map in 
Fig. |2] and listed in Table 1. Three faintest variables that do not have numbers in the 
WEBDA database were designated as VI, V2, and V3. In addition, we list in Table 
1 the four other suspects (NGC 1502-61, 37, 44 and 5) that appeared to be constant 
in our photometry. For these stars, the detection threshold in mmag is given in the 
fourth column of Table 2 instead of the variability ranges. 
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Table 1 

The variable stars and the suspected variable stars in NGC 1502 



WEBDA 


V 


V-I C 


AV, AI C 


Period 


Cluster 


Type of variability, 


number 


[mag] 


[mag] 


[mag] 


[d] 


member 


notes 


1 


6.91 


0.65 


0.08,0.07 




yes 


unknown (NSV 1458) 


2 


6.91 


0.61 


0.30,0.28 


2.698415(7) 


yes 


EA (SZ Cam) 


26 


9.63 


0.74 


0.007,0.006 


0.09611475(8) 


yes 


P Cep (NSV 15899) 


23 


10.69 


0.81 


0.08,0.09 




yes 


Be 


43 


11.36 


0.72 


0.016,0.014 


0.68613(9) 


yes 


Ell or SPB 


7 


11.49 


0.76 


0.16,0.16 


1.984476(6) 


yes 


EA 


27 


12.42 


0.81 


0.07,0.08 


3.3085(4) 


yes 


a 2 CVn 


89 


14.65 


1.47 


0.05,0.03 


0.74266(11) 


uncert. 


unknown 


VI 


15.45 


1.31 


0.52,0.44 


0.2992291(4) 


no 


RRc 


V2 


15.47 


1.61 


0.11,0.06 


1.00695(12) 


uncert. 


XEri? 


V3 


15.91 


1.42 


— ,0.28 


1.76845(13) 


uncert. 


EW? 


61 


9.56 


0.63 


<2.7,<3.3 




yes 


est (NSV 15901) 


37 


9.65 


0.70 


<1.5,<1.6 




yes 


est (NSV 15898) 


44 


10.76 


0.75 


<2.4,<2.6 




yes 


est (NSV 1456) 


5 


11.41 


0.72 


<6,<6 




yes 


est (NSV 1457) 



4.1. p Cep star NGC 1502-26 (NSV 15899) 

As noted above, NGC 1502-26 was the best (3 Cep candidate in NGC 1502, dis- 
covered by Delgado et al. (1992) and later confirmed by Steslicki (2006) using data 
obtained with the old camera. The observations carried out with the new camera 
lowered the detection threshold, but the mode with frequency of 10.404 d _1 re- 
mained the only one detected in the photometry of this star. It was found in the B , 
V , and Iq -filter data, but not in U ; apparently, the accuracy of the U photometry 
was insufficient to allow detection. The Fourier amplitude spectrum for NGC 1502- 
26 V -filter combined 2005-2008 data is shown in Fig. [3] Decreased noise in the 
low-frequency domain is due to detrending. It was done in the following way. Af- 
ter removing the periodic signal, the nightly means were calculated using residuals. 
These means were subsequently interpolated by a cubic spline and such smoothed 
curve was subtracted from the original data. 

The two radial-velocity observations obtained for this star by Huang and Gies 
(2006a) differ by about 35 km/s. They concluded therefore that the star is a single- 
lined spectroscopic binary. This is still a possibility, but a pulsation is a much more 
likely cause of the difference. The projected rotational velocity derived by the same 
authors amounts to 174 ± 8 km/s; the high rotational velocity of NGC 1502-26 is 
consistent with the MK spectral type, B 1 .5 Vn, reported by Hoag and Applequist 
(1965). 
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Fig. 3. Top: Fourier amplitude spectrum of the combined 2005-2008 V -filter data of NGC 1502-26. 
Bottom: The same after prewhitening with the mode having frequency of 10.404230 d _1 . 

The parameters of the sine-curve fit for the single mode with a frequency of 
10.404230 ± 0.000009 d" 1 , found in NGC 1502-26, are given in Table 2. The 
columns are self-explanatory; we only note that observations with the old and new 
camera are designated as 'O' and 'N', respectively. Numbers in the parentheses 
are the standard deviations with the leading zeroes suppressed. The other period- 
icity suggested by Delgado et al. (1992) for NGC 1502-26, having frequency of 
15.7 d _1 , was not found in our data. 



Table2 

Parameters of the sine-curve fit for NGC 1502-26. T = HJD 2 453 000 



Filter 


Nobs 


Semi-amplitude 


T m -dx — T(j 


S/N 


Ores 


(camera) 




[mmag] 


[d] 




[mmag] 


B (N) 


762 


3.39(28) 


1487.2497(12) 


6.2 


5.40 


V(0) 


2402 


3.38(15) 


394.9056(07) 


12.6 


5.01 


V (N) 


1498 


3.80(19) 


1360.6650(08) 


8.1 


5.15 


V(0+N) 


3900 


3.63(12) 


765.8114(05) 


14.5 


5.07 


/c(O) 


652 


3.30(38) 


361.6477(20) 


5.6 


7.16 


fcCN) 


1693 


3.09(15) 


1489.8437(08) 


8.5 


4.46 


7 C (O+N) 


2345 


3.12(16) 


1176.1251(08) 


10.0 


5.35 



4.2. Two Periodic B-Type Variables 

Two other early-type members of the cluster appear to be periodic variables. 
NGC 1502-27 was observed only with the new camera. The variability of this star 
is slightly non-sinusoidal in shape (Fig. HJ). The three 2005 nights, on which only 
V observations were made, indicate that the amplitude was smaller in 2005 than 
in 2008. The location of NGC 1502-27 in the color-color diagram (Fig. O is 
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Phase Phase 



Fig. 4. Left: The averaged differential B (open circles), V (dots) and Iq (plus signs) observations 
of the a 2 CVn-type star NGC 1502-27 phased with the period of 3.3085 d. Right: The same for the 
SPB suspect, NGC 1502-43, phased with a period of 0.68613 d. In both panels offsets were applied 
in order to separate light curves in different bands. 



consistent with a cluster late B-type star, that matches the Harvard classification 
of B8 given by Zug (1933). The star was found to be an Ap star by means of 
the Geneva photometry by North and Cramer (1981). Its peculiarity was recently 
confirmed by Paunzen et al. (2005) who used Aa photometry. This, as well as 
the period of variability («3.3 d), show that the star is an a 2 CVn-type variable. 
The star appears in the catalogue of Ap and Am stars as Renson 6525 (Renson and 
Manfroid 2009). 

The only SPB candidate we found is NGC 1502-43, a B6-type star according 
to Dombrovskii and Hagen-Thorn (1964). The light variability of this star can be 
represented by a single sinusoidal term with a period of about 0.687 d (Fig. HJ) 
and semi-amplitudes equal to 7.0 ± 0.5, 7.5 ± 0.6 and 5.9 ± 0.4 mmag in B, V , 
and Iq , respectively. Huang and Gies (2006a) report two measurements of radial 
velocity for this star; they differ by about 100 km/s. In view of the small photo- 
metric amplitudes, it seems rather unlikely that the radial velocity difference is due 
to pulsations. Although pulsation as the cause of variability cannot be rejected, the 
more plausible explanation of the observed light and radial velocity variations is 
that NGC 1502-43 is a non-eclipsing spectroscopic binary showing low-amplitude 
photometric modulations due to proximity effects. In other words, it is rather an 
ellipsoidal variable than an SPB star. If this were indeed the case, the period given 
in Table 1 would be equal to half the orbital period. 

4.3. SZ Camelopardalis 

SZ Cam (NGC 1502-2, HD 25638, BD + 61°676N) is the brightest member 
of the cluster and the B component of the visual double ADS 2984 (separation 
18")- The star was discovered as a spectroscopic binary by Plaskett (1924). Sub- 
sequently was found to be an eclipsing binary (Guthnick and Prager 1930), but 
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Fig. 5. Phase diagram of the U, B, V and Iq observations of the eclipsing system SZCam 
(NGC 1502-2). The V and Iq observations from the old CCD camera are shown with plus signs. 
Phase 0.0 corresponds to the primary minimum. For clarity, shifts between light curves in different 
bands were applied. 

the first comprehensive photometric study of the system was published by Wes- 
selink (1941). This author was the first to derive the true orbital period, equal to 
2.698 d. These early works were followed by many other studies (Kopal and Shap- 
ley 1956, Heintze and Grygar 1970, Kitamura and Yamasaki 1972, Budding 1974, 
1975, 1977, Mardirossian et al. 1980) in which the authors claimed that the system 
has a semi-detached configuration with a mass ratio of 0.2-0.3. Consequently, the 
increase of the orbital period found by Chochol (1977, 1980) was interpreted in 
terms of mass exchange between components. However, observations made with 
a higher spectral resolution (Mayer et al. 1994, Lorenz et al. 1998, Harries et al. 
1998) allowed a better separation of the lines of the components resulting in im- 
proved, much higher mass ratio, close to 0.7. The new solutions showed the system 
to be detached. Earlier, Mayer (1987) interpreted the changes of the orbital period 
of SZCam in terms of the light-time effect thus suggesting also the existence of 
a tertiary component in the system. The spectral lines of the tertiary were indeed 
found in the spectra obtained by Mayer et al. (1994) and the tertiary appeared to be 
itself a single-lined spectroscopic binary having an orbital period of 2.797 d (Mayer 
et al. 1994, Gorda 2002, Michalska et al. 2007). The wide pair orbits the center 
of mass of the quadreple system in an eccentric (e = 0.78) orbit with a period 
of about 50 yr (Mayer et al. 1994, Lorenz et al. 1998, Gorda et al. 2007, Gorda 
2008). The pair was also resolved by means of speckle interferometry (Mason et 
al. 1998, Gorda et al. 2007). A combination of the visual and spectroscopic (from 
the study of the light-time effect) orbits allowed to derive the distance to SZ Cam, 
1125 ± 135 pc (Gorda et al. 2007), in a reasonable agreement with the determina- 
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tions of cluster distance based on photometric methods and spectroscopic parallax. 
The total mass of the NGC 1502-2 quadruple system exceeds 50 M Q . The eclipsing 
pair is classified as 09 IV + B0.5 V (Lorenz et al. 1998). 

The UBVIc light curves of SZ Cam from our observations are shown in Fig. [5] 
We have derived 14 epochs of minimum light which are listed in Table 3. They 
were obtained by a least-squares fit of a smoothed light-curve to data subsets cov- 
ering an epoch of minima. The method uses the whole light curve, thus allowing 
determination of an epoch of minimum light even if it was only partially covered 
by observations. The derived epochs of minimum light were used to obtain the 
following ephemeris: 

Min I = HJD 2453355.0682(21) + 2.698415(7) x E, 

where E is the number of the orbital cycles elapsed from the initial epoch and the 
numbers in parentheses are the rms errors with leading zeroes omitted. We would 
like to point out, however, that the orbital period of SZCam varies due to the light- 
time effect, so that over a longer time scale the ephemeris given by Gorda et al. 
(2007) should be used. 

Table3 

Epochs of minimum light for SZ Cam and NGC 1502-7 



SZ Cam NGC 1502-7 



Filter 


Minimum 


Tmax -HJD 2453000.0 


Filter 


Minimum 


W -HJD 2453000.0 


V 


sec. 


356.4167(16) 


V 


prim. 


350.7594(05) 


Ic 


sec. 


356.4090(41) 


V 


prim. 


408.3038(05) 


V 


sec. 


410.3832(09) 


Ic 


prim. 


408.3025(14) 


Ic 


sec. 


410.3785(14) 


V 


prim. 


410.2904(05) 


V 


sec. 


429.2812(17) 


Ic 


prim. 


410.2907(11) 


k: 


sec. 


429.2747(64) 


V 


prim. 


622.6299(07) 


V 


prim. 


649.2037(11) 


B 


prim. 


1521.6051(15) 


u 


prim. 


1469.5114(33) 


V 


prim. 


1521.6094(16) 


B 


prim. 


1469.5135(08) 


Ic 


prim. 


1521.6089(07) 


V 


prim. 


1469.5115(09) 


B 


prim. 


1553.3365(35) 


Ic 


prim. 


1469.5118(12) 


V 


prim. 


1553.3365(19) 


U 


sec. 


1511.3323(45) 


Ic 


prim. 


1553.3462(10) 


B 


sec. 


1511.3457(20) 


B 


prim. 


1555.3215(14) 


Ic 


sec. 


1511.3423(20) 


V 


prim. 


1555.3219(15) 








Ic 


prim. 


1555.3249(09) 



4.4. Other Eclipsing Binaries 

NGC 1502-7 appear to be another EA-type eclipsing binary. It is also a likely 
member of the cluster. The times of the primary minimum are listed in Table 3. 
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0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 
Orbital phase 



Fig. 6. The BVIq (from top to bottom) observations of NGC 1502-7 phased with the period of 
1.984476 d. Some offsets were applied to separate light curves in different bands. Phase 0.0 corre- 
sponds to the primary minimum. 

They were derived in the same way as for SZ Cam and resulted in the following 
ephemeris: 

Min I = HJD 2453350.7565(20) + 1.984476(6) x E. 

The light curves of NGC 1502-7 (Fig. [6]) show a 0.1 8 mag deep primary minimum 
with a flat bottom and a very shallow («0.02 mag) secondary minimum. This 
shows that the secondary component is much fainter and smaller than the primary. 
The system is detached; the components are probably main-sequence stars. Un- 
fortunately, the two radial- velocity measurements for NGC 1502-7 given by Huang 
and Gies (2006a) are separated by almost exactly one orbital period; this does not 
allow to estimate the range of radial velocities for the system's primary. 



5.2 - 




0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 

Phase 

Fig. 7. Phase diagram of the 2008 1q observations of V3, possible eclipsing EW-type system. Phase 
0.0 was chosen arbitrarily. 

Another variable star that might be an eclipsing binary is V3, a relatively faint 
star located close to the cluster center. The proximity of the brightest stars in the 
field affected the photometry of V3 so much that the variability can be clearly seen 
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only in the 2008 Iq data and barely in the 2008 V -filter photometry. The shape 
of the light curve (Fig. resembles W UMa-type variability. There is also a weak 
evidence for the O'Connell effect. However, the period (1.77 d) is rather long for 
this type of eclipsing binary. Therefore, we cannot exclude the possibility that the 
true period is half that given in Table 1 and that the cause of variability is other than 
eclipses. 




0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 

Phase 



Fig. 8. The averaged differential B (open circles), V (dots) and Iq (plus signs) observations of the 
background RR Lyr star VI phased with the period of 0.2992291 d. Offsets were applied to separate 
light curves in different bands. 

4. 5. The RR Lyrae Star VI 

VI is clearly an RRc-type star judging from the period, shape of the light curve 
and increasing amplitude when going from Iq to B (Fig. [8). The star is definitely 
not a member because the cluster is much too young to contain such an evolved star. 
Adopting the instrinsic (V —Ic)o = 0.35 mag for an RRc star (see e.g., Kopacki 
2001), we getE(V — /c) ~ 0.99 mag corresponding to E(B—V) ~ 0.8 mag if color 
excess ratio given by Winkler (1997) is applied. Adopting Ry = 2.6 (see Section 
2) we get Ay « 2.1 mag. The absolute magnitude of the star cannot be estimated 
accurately because its metallicity is not known. However, assuming a typical value 
of [Fe/H] = — 1.5 and using the relation given by Olech et al. (2003), we get an 
estimate of the true distance modulus, 12.8 mag, and a distance of 3.5 kpc. VI 
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is therefore a background object, only slightly more reddened than the cluster, but 
located at a much larger distance than NGC 1502. 



<d 10.5 

T3 

Z3 
H — ; 

O) 10.6 

CO 



_p 10.7 

> 
CO 

z> 



10.8 



NGC 1502-23 



o v C 



800 900 1000 1100 2000 
HJD 2453000+ 

Fig. 9. UBVIc (from top to bottom) light curves of NGC 1502-23. The following offsets in magni- 
tudes have been applied: U : —0.78 mag, B : —0.75 mag, V : no offset, Iq '■ +0.9 mag. Each point 
corresponds to a nightly mean. 



4.6. Variable Stars Showing H(X Emission 

Fast rotation and the presence of circumstellar disks in Be stars result in vari- 
ability on different time scales, both regular and irregular. Of the stars with emis- 
sion we detected by measuring the a index (Section 6), two appeared to be variable. 
NGC 1502-23 showed irregular behavior (Fig. [9]> typical of Be stars. Its MK spec- 
tral type as determined by Hoag and Applequist (1965) is B3 V; the location in the 
two-color diagram (Fig. [[4]) is consistent with this classification. 



CO 

E, 4.3 

CD 

■D 4.4 
o> 4.5 

CO 

I 4.6 

CO 

| 4.7 

CD 

£ 4.8 
Q 



0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 

Phase 

Fig. 10. VIq light curves of V2 phased with the period of 1.00695 d. 

Another star with Ha emission, V2, is also a variable. The dominant periodic- 
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ity has a period of about 0.503 d or, more likely, twice this value, 1.007 d (Fig.fTOl. 
The proximity of the period to 1 d makes difficult the derivation of an unambiguous 
period from the combined 2005 and 2008 data. In addition to the period given in 
Table 1, another variation was found in the Iq -filter data, with a period of 1.966 d, 
i.e., about twice as long as the dominant period. Fig. [10] shows the Vic data for 
this star, averaged over 0.0025-day intervals and phased with the 1.007-day period. 
As the spectral type of this star is not known, it is difficult to assign the type of 
variability, but the period and the Ha emission indicate a X Eri type. 

4. 7. Other Variables 

As noted in Section 2, NGC 1502-1 = ADS 2948A = NSV 1458 was claimed to 
be variable by Guthnick and Prager (1930). The star has a spectral type of B0ITIII 
(Morgan et al. 1955) and is known to have variable radial velocity (Plaskett 1924, 
Lorenz et al. 1998, Hohle et al. 2009). Unfortunately, the photometric variability of 
NGC 1502-1 (Fig. ITTb is rather erratic. We also analyzed archival radial velocities 
in the hope of finding periodic variations. Unfortunately, due to the scarcity of the 
radial velocity data, the results were inconclusive. 



w 

CD 
T3 

Z3 
■ 4— ; 

'c 

D) 
03 

E 

_o 

> 

CO 
=> 




800 



900 



2000 



1000 1100 
HJD 2453000+ 

Fig. 11. UBVIq (from top to bottom) light curves of NGC 1502-1. The following offsets in magni- 
tudes were applied: U : —0.25 mag, B : —0.55 mag, V: no offset, Iq '■ +0.72 mag. 



Finally, periodic low-amplitude variations were found in the photometry of 
NGC 1502-89 (Fig. [12]>. The period is equal to 0.74266 d. However, the spec- 
tral type of this star is not known and consequently no final type of variability can 
be assigned. 

As mentioned in Section 2, four other stars in the field, NGC 1502-61, 37, 44 
and 5, were claimed to be photometrically variable. We do not find any evidence 
for their variability with semi-amplitudes higher than the values given in the fourth 
column of Table 1. Unlike for the variable stars, the detection thresholds for these 
stars, corresponding to S/N = 4 in the Fourier periodogram, are given in mmag. 
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Q ■ Al c + 1.25 



0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 

Phase 

Fig. 12. The differential 2008 V (dots) and Iq (plus signs) observations of NGC 1502-89 phased 
with the period of 0.74266 d. 

The time-series photometry of all variable stars we found in NGC 1502 has 
been deposited in the Acta Astronomica Archive. 



As pointed out in Section 1, the UBV photomeuy of NGC 1502 was provided 
by many authors, but except for Hohle et al. (2009), no Iq photometry has been 
made so far. In order to obtain standard (V — Iq) colors, we carried out Vic ob- 
servations of NGC 1502 together with the Landolt (1992) field PG 0918 +029 on 
a single clear night of April 20, 2005. From these data, we derived the (V — Iq) 
colors of 29 stars in the central part of the field. These stars were subsequently used 
as secondary standards for the transformation of the instrumental (V — Iq) colors 
obtained with the new camera. The transformation equation was the following: 



where a denotes standard deviation of the fit while numbers in parentheses, the 
rms errors of the coefficients with leading zeroes omitted. In the transformation 
equations we use lowercase letters for the instrumental and uppercase for the stan- 
dard magnitudes. The V magnitudes, as well as (U — B) and (B — V) colors, were 
obtained by means of a transformation to the standard system using the UBV obser- 
vations of Hoag et al. (1961), both photoelectric and photographic. The following 
transformation equations were derived: 



5. UBVIq Photometry 



(V-I c ) =0.745(17) x (v-ic) +0.739(06), o 



0.029 mag 



V = v - 0.134(18) x (v - i c ) + 9.628(07), a 



0.043 mag, 



(U-B) = 1.049(48) x (u — — 0.229(26), a 
(B-V) = 1.350(30) x (b-v) + 0.599(08), a 



0.048 mag. 



0.122 mag 
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On average, about 50 stars were used to derive each transformation. The coeffi- 
cients for color terms considerably different from 1 for transformations to (B — V) 
and (V — 7c) and different from for transformation to V in the above equations 
are the consequence of the fact that our instrumental v band is shifted towards 
shorter wavelengths with respect to the standard V . In total, V magnitudes and 
(V —7c) colors were derived for 191 stars in the field, (B — V) colors for 119 stars 
and (U — B) colors, for 55 stars. The resulting magnitudes and colors were used 
to plot the stars from NGC 1502 in the color-magnitude and color-color diagrams 
(Figs. [13] and IT4T >. The photometry is given in Table 4. In addition, astrometric 
transformation of mean stellar positions derived from images was made in order to 
obtain equatorial coordinates of the observed stars. As a reference catalogue, we 
used the USNOB1.0 catalogue (Monet et al. 2003). A standard deviation of the 
transformation amounts to about 0.1". The positions of stars are also given in Ta- 
ble 4. However, we provide here only the entries for variable stars. The full version 
of Table 4 is available in electronic form from the Acta Astronomica Archive. 



7 
8 
9 

10 
11 

O) 
CO 

£ 12 
> 

13 
14 
15 
16 



(m - M) = 10.0 mag 
E(V- l c ) = 0.90 mag 
= 1 .82 mag 




0.5 



1.0 1.5 
(V-I c ) [mag] 



2.0 



Fig. 13. The V vs. (V — Iq) color-magnitude diagram for NGC 1502. The eleven variable stars are 
shown as encircled dots and labeled. The three isochrones were taken from the models of Bertelli et 
al. (1994) for Z = 0.019. They correspond to log(age/yr) = 6.6, 7.0 and 7.3, i.e. 4, 10 and 20 Myr. 
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0.0 0.2 0.4 0.6 0.8 1.0 1.2 
(B - V) [mag] 

Fig. 14. The (U — B) vs. (B — V) color-color diagram for the brightest stars in the observed field 
of NGC 1502. The dashed line shows the intrinsic color-color relation for main- sequence stars as 
given by Caldwell et al. (1993). The same relation for reddened stars with E(B — V) = 0.60, 0.75 
and 0.82 mag is plotted with continuous lines. It was assumed that E(U — B) /E(B — V) = 0.724 + 
0.021 xE(B-V) (Turner 1989). 



Table4 

UBVIq and Ha photometry and coordinates of stars in the field of NGC 1502 



Star 


R. Asc. (2000.0) 


Decl. (2000.0) 


V 


(V-fc) 


(B-V) 


(U-B) 


a 


1 


4 h 07 m 5 1.38 s 


+62°19'48.4" 


6.910 


+0.653 


+0.490 


-0.373 


2.113 


2 


4 h 07 m 49.27 s 


+62°19'58.7" 


6.915 


+0.611 


+0.445 


-0.395 


2.136 


26 


4 h 07 m 44.08 s 


+62°18'04.2" 


9.630 


+0.744 


+0.590 


-0.232 


2.091 


23 


4 h 07 m 5 1.75 s 


+62° 16' 29.7" 


10.692 


+0.814 


+0.654 


-0.058 


2.056 


43 


4 h 07 m 49.75 s 


+62°22'25.2" 


11.357 


+0.716 


+0.576 


-0.085 


2.193 


7 


4 h 07 m 54.88 s 


+62°20'17.7" 


11.488 


+0.762 


+0.627 


+0.031 


2.232 


27 


4 h 07 m 40.81 s 


+62°15'40.7" 


12.420 


+0.814 


+0.579 


-0.017 




89 


4 h 08 m 08.98 s 


+62°15'36.0" 


14.651 


+ 1.474 


+ 1.453 






VI 


4 h 07 m 23.10 s 


+62°24'14.3" 


15.450 


+ 1.311 


+ 1.176 






V2 


4 h 07 m 59.92 s 


+62° 17' 17.6" 


15.470 


+ 1.607 


+ 1.617 




1.898 


V3 


4 h 07 m 50.82 s 


+62° 19' 08.2" 


15.912 


+ 1.419 









As can be seen from Fig.[T3] the cluster main sequence is well defined, although 
some scatter (f»0.2 mag) in (V — Ic) due to differential extinction can be seen. 
In order to estimate the age of the cluster, we adopted the distance of 1 kpc, i.e., 
(m—M)o = 10.0 mag and the mean color excess E(B — V) of 0.70 mag, translating 
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into E(V — Iq) ~ 0.90 mag and a total absorption in V , Ay = 1.82 mag (Ry = 2.7 
was assumed, see Section 2). All three isochrones we show in Fig. |T3]fit the cluster 
main sequence equally well and it is the position of the two brightest stars that 
allows us to constrain the age as the turnoff point of the cluster is poorly populated. 
The comparison allows us to estimate the age of the cluster for 10 ± 5 Myr, in a 
full agreement with previous determinations. 



6. Ha Photometry 

As for most clusters surveyed within our program, we made Ha photometry 
of NGC 1502 using two filters centered at Ha, narrow and wide. The resulting 
reddening-free a index can be used to find stars with emission at Ha, including 
Be stars. The Ha observations for NGC 1502 were carried out with the old camera 
and therefore the a index could be derived only for the brightest stars in the central 
part of the cluster. The indices are shown in Fig. [15] and given in the last column of 
Table 4. 



1.9 

2.0 
S 3 2.1 
— 2.2 

a 

2.3 
2.4 
2.5 

7 8 9 10 11 12 13 14 15 
V [mag] 

Fig. 15. The a index plotted against V magnitude for the brightest stars in the field covered by the 
observations made with the old CCD camera. The four crosses denote relatively faint stars located in 
the vicinity of the two brightest stars in the field, NGC 1502-1 and 2. Open circles are stars that could 
be identified as non-members from the color-magnitude diagram (Fig. [Oi l. The solid line shows 
calculated location of luminosity class V stars that are cluster members; see text for explanation. 




Kolaczkowski et al. (2004) showed how a changes for stars of different spec- 
tral types and luminosity classes. We adopted these calculations, but expressed 
the calculated values of a as a function of absolute magnitude using the absolute 
magnitude vs. spectral type calibration of Schmidt-Kaler (1982). Then, adopting 
the apparent distance modulus of NGC 1502, (m — M) = 11.8 mag, we plotted the 
resulting a vs. V relation in Fig. [TBI (continuous line). The relation has a maxi- 
mum (note the reversed ordinate scale) for stars with early A spectral type because 
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these stars have the strongest Ha lines. The points follow the relation because the 
majority of bright stars belongs to the cluster and have no emission at Ha. Stars 
located well below the relation, like NGC 1502-11, are non-members; for those 
located above it, there are two possibilities: either they are non-members or have 
emission at Ha. The two possibilities are not mutually exclusive, i.e., we can have 
a non-member with Ha in emission. However, all stars having a < 2.0 mag could 
be regarded as stars with Ha emission. We have three such stars, including the 
variable star V2. Some non-members can be identified from the color-magnitude 
diagram (Fig.fTBl where they are located clearly off the cluster main sequence. Four 
such stars were plotted in Fig. [15] as open circles. However, there are still several 
stars that might have weak emission in Ha. These are: NGC 1502-26, the [3 Cep 
star (see Section 4.1), NGC 1502-29, 25, 51 and 88. 



7. Conclusions 



The open cluster NGC 1502 contains only a single low-amplitude p Cep star 
(NGC 1502-26). This appears to be in a contrast to NGC 6910 (Kolaczkowski et 
al. 2004, Pigulski et al. 2007) which is not very rich in stars either but harbors as 
many as seven [3 Cep stars. Whether the difference can be interpreted in terms 
of a difference in metallicity is not known. Unfortunately, the only metallicity 
determination that can be related to NGC 1502 is that from Stromgren photometry 
by Eggen (1985) for HD 25056 in Cam OBI. He reported [Fe/H] = 0.16 dex for 
this star. The value of [Fe/H] is rather uncertain. In addition, the star might not be 
related to NGC 1502 at all. Definitely, a spectroscopic determination of abundances 
of stars in the cluster are highly desirable. 

SZ Cam remains one of the most interesting stars of the cluster being an ex- 
cellent case of a massive hierarchical system that are found frequently in the cores 
of open clusters. Thorough studies of such systems may shed light on the role of 
binaries in massive star formation in general. SZ Cam was also succesfully used to 
derive the distance to the cluster (Gorda et al. 2007) and we may hope that further 
observations will improve this determination. The relatively bright EA-type system 
(NGC 1502-7) we discovered might be also used for this purpose. 
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